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ABSTRACT 

The  development  of  electromechanical  actuators  (EMAs)  is  the  key  technology  to  build  an  all¬ 
electric  aircraft.  One  of  the  greatest  hurdles  to  replacing  all  hydraulic  actuators  on  an  aircraft 
with  EMAs  is  the  acquisition,  transport  and  rejection  of  waste  heat  generated  within  the  EMAs. 
The  absence  of  hydraulic  fluids  removes  an  attractive  and  effective  means  of  acquiring  and 
transporting  the  heat.  To  address  thennal  management  under  limited  cooling  options,  accurate 
spatial  and  temporal  information  on  heat  generation  must  be  obtained  and  carefully  monitored. 

In  military  aircraft,  the  heat  loads  of  EMAs  are  highly  transient  and  localized.  Consequently,  a 
FEA-based  thermal  model  should  have  high  spatial  and  temporal  resolution.  This  requires 
tremendous  calculation  resources  if  a  whole  flight  mission  simulation  is  needed.  A  lumped  node 
thennal  network  is  therefore  needed  which  can  conectly  identify  the  hot  spot  locations  and  can 
perform  the  calculations  in  a  much  shorter  time.  The  challenge  in  forming  an  accurate  lumped 
node  thermal  network  is  to  determine  all  the  suitable  thennal  resistances  and  capacitances  of  the 
thennal  network. 
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In  this  paper  we  present  an  FEA-based  lumped  node  network  and  its  simulation  of  a  mission 
profile.  This  model  is  based  on  a  detailed  FEA  model  to  locate  the  hot  spots,  to  detennine  the 
network  parameters  and  to  verify  its  effectiveness.  The  model  can  also  deal  with  the  nonlinear 
behavior  of  the  EMA  system  introduced  by  phase-change  materials  (PCM)  if  thermal  energy 
storage  is  needed,  and  temperature-dependent  magnetic  properties.  This  model  can  also  be 
incorporated  into  lumped  node  magnetic  and  electric  model  to  develop  a  full  multi-physics, 
multi-scale  simulation  engine.  This  engine  can  accurately  analyze  the  complete  EMA  system  in  a 
systematic  scale  and  whole-mission  duration. 

INTRODUCTION 

The  application  of  EMA  on  military  aircraft  poses  a  challenge  in  the  design  of  the  whole  driving 
train,  from  the  motor-controller  to  the  gearbox.  To  successfully  replace  the  hydraulic 
counterparts,  the  EMAs  have  to  provide  reliable  driving  force  in  whole  mission  profile  under  the 
thennal,  power,  mass  and  size  constraints.  In  addition,  the  EMA  performance,  reliability  and 
power  consumption  depend  on  the  temperatures  of  the  key  components,  including  the  power 
electronics  and  electric  motor.  Waste  heat,  if  not  properly  managed,  could  cause  high 
temperatures  of  the  motor-driver,  motor  and  gearbox,  leading  to  perfonnance  degradation,  and 
creating  even  more  power  losses,  and  hence  would  increase  the  temperatures  even  more.  A 
complete  multi-physics  simulation  study  for  a  realistic  mission  profile  is  essential  to  fully 
understand  the  perfonnance  of  the  EMAs  on  aircraft. 

Finite  Element  Analysis  (FEA)  is  widely  used  in  today’s  mechanical,  electronic  and 
aerodynamic  designs.  However,  direct  application  of  FEA  to  analyze  the  whole  EMA  system 
would  require  tremendous  amount  of  computational  power.  An  estimate  of  one-hour  duration  of 
FEA  transient  thermal  simulation  of  an  electric  motor  requires  two  terabytes  of  hard  drive  disk 
space  and  720  hours  of  CPU  time  on  a  typical  quad-core  computer.  The  whole  EMA  model  is 
more  complex  than  the  electric  motor  and  consequently  requires  even  more  computational 
resources.  A  multi-physics  FEA  simulation,  including  electronics,  magnetics  and  stress 
calculations  would  make  such  a  simulation  impractical,  especially  when  real-time  information  is 
needed. 

Another  approach  is  to  develop  a  lumped  node  thermal  network  to  represent  the  temperature  of 
every  solid  part  of  the  EMA.  The  thermal  resistances  and  capacitances  between  the  nodes  can  be 
treated  as  electrical  resistances  and  capacitances.  Hence  the  temperature  of  every  node  can  be 
solved  as  the  voltage  in  this  equivalent  network  (Table  1).  This  approach  is  well  developed  in 
motor  design  industry  [1,2].  Some  commercial  motor  design  software  package  has  already 
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included  the  thermal  network  simulation,  i.e.  Motor-CAD  [3],  Some  studies  have  been  made 
with  FEA  analysis  and  experimental  testing  have  shown  that  such  an  approach  is  valid  [4]. 
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Equations  1  and  2  show  how  to  calculate  the  R  and  C  values  for  simple  one-dimensional 
geometry.  Accurate  estimates  for  R  and  C  values  are  not  possible  for  complicated  geometries 
such  as  an  electric  motor.  The  traditional  lumped  node  thermal  network  is  based  on  a  “forward” 
direction  modeling  process.  In  this  process,  with  all  the  material,  geometry  and  construction 
information  provided,  the  thennal  resistance  R0  and  capacitance  C9  are  calculated  based  on 
empirical  and  theoretical  relations.  For  example,  when  the  winding  wire  type,  winding  structure, 
potting  material,  slot  width  and  teeth  thickness  are  given,  the  thermal  resistance  from  winding  to 
stator  and  from  winding  to  rotor  can  be  calculated.  Because  this  type  of  approach  requires  large 
number  of  empirical  relations  to  achieve  high  accuracy,  a  “reverse”  modeling  process  is 
introduced  in  this  paper.  A  detailed  3-D  solid  model  of  a  target  motor  is  constructed  and 
imported  to  an  FEA  software  like  ANSYS  [5]  to  perform  steady-state  and  transient  simulation. 
With  these  results  we  can  estimate  and  select  the  values  for  the  thermal  resistances  and 
capacitances.  The  advantage  of  this  method  is  that  we  can  evaluate  the  fidelity  of  the  lumped 
node  model  with  a  real  motor,  add  or  reduce  nodes  to  increase  the  model  accuracy  and 
efficiency,  and  estimate  the  maximum  error  between  the  node  temperature  and  maximum 
temperature  in  real  motor. 

This  procedure  can  also  be  extended  to  model  the  gear-box,  motor-driver  and  drive-train,  and 
even  include  the  aircraft  wing  surfaces  and  frames.  The  thennal  network  can  also  be 


3 


incorporated  into  a  multi-physics  model  to  simulate  the  electrical,  thermal  and  mechanical 
performance  of  the  whole  EMA  and  its  supporting  structure.  Once  the  proper  thermal  resistances 
and  capacitances  are  selected,  this  simulation  engine  can  be  used  with  various  time  dependent 
boundary  conditions  during  the  whole  mission  duration,  including  the  air  temperature,  aircraft 
speed,  altitude  and  sunlight.  The  computational  requirement  of  such  a  simulation  is  negligible 
compared  to  the  FEA  simulation. 

MODELING 

Figure  1  is  the  3-D  model  of  a  typical  Permanent  Magnet  Synchronous  Machine  (PMSM)  servo 
motor.  This  design  features  a  12-slot  stator  and  a  10-pole  rotor  (Figure  2).  This  motor  is  designed 
to  input  10  hp  electrical  power,  with  a  efficiency  of  91.8%.  The  total  power  loss  in  the  motor  is 
61 1W.  In  the  computation  domain  (Figure  1)  the  power  loss  is  153W.  The  distribution  of  the 
power  loss  in  the  motor  components  is  listed  in  Table  3.  In  Figure  1  the  parts  are  numbered  in 
the  same  way  as  in  Table  2.  These  numbers  are  also  the  same  as  those  in  the  lumped  node 
network  shown  in  Figures  6  and  7. 


Figure  1.  3-D  model  of  the  PMSM  motor  design  and  computation  domain. 
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Figure  2.  Cross-section  of  motor  and  rotor. 

Due  to  the  symmetry,  the  geometry  of  the  model  is  divided  into  a  quarter  section  and  imported 
into  ANSYS  for  the  FEA  simulation.  A  total  thermal  load  given  in  Figure  4  and  with  heat 
distribution  shown  in  Table  3  are  assumed.  An  external  forced  convective  heat  transfer 
coefficient  of  100  W/m2.K  applied  to  the  case  of  the  motor.  Figures  3  and  5  show  the  steady- 
state  temperature  field  and  the  transient  temperature  response  inside  the  windings  and  rotor 
magnet. 
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Figure  3.  Steady  state  thermal  simulation  result  (forced  convection). 
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Figure  4.  Step  heat  load  of  the  motor 


Figure  5.  Transient  temperature  response  of  the  load  in  Figure  4 
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ANSYS  (or  any  other  FEA  software)  generally  gives  temperature  infonnation  in  the  whole 
calculation  domain.  However,  in  practice,  only  several  locations  or  components  are  needed  to 
assess  the  thermal  behavior  of  the  motor.  For  the  problem  at  hand,  we  focus  only  on  the 
temperatures  of  the  winding,  stator,  magnet  and  bearing.  It  is  noted  that  for  every  component,  the 
temperature  is  not  unifonn.  The  temperature  range  can  be  characterized  by  the  maximum  and 
minimum  temperatures  within  the  particular  component.  For  the  copper  winding,  these  two 
temperatures  are  quite  close  because  of  the  excellent  thermal  conductivity  of  copper.  While  in 
the  stator  these  two  temperatures  differ  by  a  few  degrees.  In  the  lumped  node  network  we  only 
use  a  single  value  to  represent  each  component,  which  generally  is  the  average  temperature.  The 
deviation  between  the  average  temperature  and  the  maximum  temperature  in  the  components 
should  be  kept  in  mind  in  the  lumped  node  simulation. 

The  node  network  is  constructed  as  shown  in  Figure  6.  The  electric  motor  is  assumed  to  have  3 
major  thennal  paths  with  other  parts  of  EMA  and  environment.  The  first  one  is  the  thermal 
resistance  R5,  which  represents  the  convection  thennal  resistance  between  motor  surface  and  the 
air  around  it.  Considering  that  the  motor  will  be  installed  in  a  small  space  inside  the  wing 
structure,  this  convection  is  likely  to  be  natural  convection  or  weak  forced  convection.  In 
ANSYS  we  assume  the  heat  transfer  coefficient  for  the  forced  convection  to  be  h=100  W/m2.K, 
which  is  a  typical  value  for  forced  convection  in  air  (Figure  3).  The  value  of  R5  can  be  easily 
changed  to  accommodate  the  actual  flow  situation  in  the  motor’s  final  installation  and  flow 
conditions. 

The  second  thennal  path  is  the  conductive  heat  transfer  from  the  motor  front  cover  to  gear-box, 
which  is  represented  as  R17.  The  third  one  is  the  conductive  heat  transfer  from  the  motor  rotor  to 
the  gear  mechanism.  The  third  one,  R7a,  could  be  a  minor  heat  flow  path  because  the  rotor  and 
bearing  loss  is  relatively  small.  But  there  is  a  possibility  that  the  gear  mechanism  has  higher 
temperature  than  the  rotor  and  conduct  heat  back  to  motor.  There  could  be  additional  thennal 
paths,  such  as  from  the  motor  back  plate  to  the  aircraft  frame,  which  are  not  shown  in  Figure  6. 
However,  incorporating  additional  thermal  paths  into  the  network  and  changing  the  lumped-node 
model  accordingly  is  simple  with  the  method  presented  in  this  paper. 

The  following  procedure  is  used  to  determine  the  R  and  C  values  in  the  lumped  node  network  of 
motor  assembly  shown  in  Figure  6.  Table  2  lists  all  the  nodal  temperatures  from  ANSYS  steady 
state  simulation  for  motor  with  natural  and  forced  convection  boundary  conditions.  In  this  Table, 
the  temperatures  with  the  presence  of  radiation  heat  transfer  is  also  considered  and  listed.  It 
should  be  noted  that  in  a  situation  where  air  cooling  is  not  effective  (such  as  natural  convection 
with  a  low  h),  the  motor  temperature  could  be  very  high  and  radiative  heat  transfer  becomes  the 
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dominant  mode  of  heat  transfer.  However  when  there  is  more  effective  convective  heat  transfer, 
the  radiation  becomes  a  minor  heat  path.  For  forced  convection,  the  Nusselt  number  for  a 
cylinder  is  detennined  from[6], 


Nud  =  CRe^Pr1/3  (3) 

where  C  and  m  are  based  on  the  Reynold’s  number.  However,  it  is  difficult  to  detennine  the 
flow  condition  and  the  characteristic  length  of  the  motor  and  drive  bay.  In  this  paper,  a  high  h 
value  of  100  to  represent  force  convection  flow  is  used  and  radiation  is  not  included.  A 
comparative  low  convective  transfers  coefficient  (h=10)  case  is  also  listed  in  table  2  to  show  the 
dominance  of  radiative  heat  transfer  in  that  situation.  The  radiation  thermal  path  can  be  easily 
included  by  adding  a  radiative  thennal  resistance  to  the  lumped  node  model  when  it  is  necessary. 
The  h  value  can  also  be  easily  altered  based  on  the  working  condition  of  the  motor  installed 
without  any  change  of  other  parameters  of  the  motor  thermal  model. 

Table  2  and  the  power  loss  table  (Table  3)  are  combined  to  solve  the  steady  state  lumped  node 
network.  (Figure  7) 
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8=0.85) 

Copper 

winding 

Copper 

1 

162 

153 

1229 

396 

Epoxy 

Epoxy 

2 

161 

152 

1228 

395 

Iron  stator 

Iron 

3 

161 

152 

1228 

395 

Case 

Aluminum 

alloy 

4 

152 

145 

1220 

385 

Case 

surface 

5 

150 

143 

1218 

382 

Front  cover 

Aluminum 

alloy 

14 

141 

131 

1210 

377 

Shaft 

Stainless 

Steel 

6 

149 

139 

1216 

384 

8 


Shaft  front 
end 

Stainless 

Steel 

7 

145 

138 

1214 

383 

Shaft  back 
end 

Stainless 

Steel 

8 

148 

138 

1214 

383 

Magnet 

NdFeB 

9 

149 

140 

1217 

385 

Airgap 

Air 

10 

160 

151 

1227 

395 

Front 

bearing 

Steel 

12 

143 

135 

1212 

381 

Rear 

bearing 

Steel 

13 

147 

137 

1212 

381 

Table  2.  Steady  state  simulation  results  with  ANSYS 


Power  loss  (W) 

Percentage  of  total 
power  loss 

IC  (copper  loss) 

144 

94.33% 

IS  (stator  loss) 

6.06 

3.97% 

IM  (magnet  loss) 

0.86 

0.56% 

IW  (windage  loss) 

0.35 

0.23% 

IBF  (front  bearing  loss) 

0.84 

0.55 

IBB  (rear  bearing  loss) 

0.54 

0.35% 

Table  3.  Component  heat  generation 
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Detailed  Lumped  node  model  Of  EMA  (transient) 


1.  copper 

2.  epoxy 

3.  stator 

4.  case  and  back  cover 

5.  case  surface 

6.  shaft  body 

7.  shaft  front  end 

8.  shaft  back  end 

9.  magnet 

10.  air  gap 

11.  connection  gear 

12.  front  bearing 

13.  back  bearing 

14.  front  cover 

15.  moving  mechanism 

16.  wing  surface 

17.  gear  box  case 


Air  flow  outside 


Motor  assembly 


Figure  6.  Lumped  node  model  of  EMA  (transient) 
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Detailed  Lumped  node  model  of  Motor  (static) 


R5 

AA/V- 


©IM(t)  R. 

Magnet 

R6  loss  pg 

6  LAA/V—* — VW 


w  Bearing 


Shaft 
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10 

Air  gap 


-A 

IW(t) 


IG(t) 


1.  copper 

2.  epoxy 

3.  stator 

4.  case  and  back  cover 

5.  case  surface 

6.  shaft  body 

7.  shaft  front  end 

8.  shaft  back  end 

9.  magnet 

10.  air  gap 

11.  connection  gear 

12.  front  bearing 

13.  back  bearing 

14.  front  cover 

15.  moving  mechanism 

16.  wing  surface 

17.  gear  box  case 


13" - WV 


8 

Shaft 
Back  end 


-WV 


Motor  assembly 


Figure  7.  Lumped  node  model  of  EMA  (steady  state) 

For  a  pure  resistance  network  with  all  the  voltages  and  currents  sources  known,  the  standard 
procedure  is  to  apply  Kirchhoff  s  Current  Law  (KCL)  equations  to  every  node.  However,  the 
number  of  nodes  is  generally  less  than  the  number  of  resistances,  which  means  there  are  not 
enough  equations  to  solve  for  the  resistances.  One  solution  is  to  add  additional  heat  flux 
conditions  from  ANSYS.  In  this  case,  heat  flux  through  RIO,  R7  and  R8  are  added.  After  that, 
the  thermal  capacitance  values  can  be  calculated  from  Equation  2.  Table  4  lists  the  R  and  C 
values  used  in  the  lumped  node  model. 
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R1 

°c/w 

0.001155 

R4 

°c/w 

0.001891 

R2 

°c/w 

0.00042 

R4a 

°c/w 

0.53 

R3 

°c/w 

0.0527 

R6 

°c/w 

0.0857 

R5 

°c/w 

0.72 

R6b 

°c/w 

1.022 

R6a 

°c/w 

1.03 

R7 

°c/w 

0.51 

RIO 

°c/w 

8.49 

R8 

°c/w 

0.47 

R12 

°c/w 

1.61 

R9 

°c/w 

20.933 

R13 

°c/w 

2.456 

R14a 

°c/w 

3.75 

Cl 

J/°C 

35.68 

C4 

J/°C 

235.1 

C2 

J/°C 

16.47 

C7 

J/°C 

14.99 

C3 

J/°C 

59.3 

C8 

J/°C 

7.08 

C6 

J/°C 

56.03 

CIO 

J/°C 

0.00176 

C9 

J/°C 

6.72 

C14 

J/°C 

62.56 

C13 

J/°C 

6.72 

C12 

J/°C 

13.18 

Table  4.  The  R  and  C  values  of  lumped  node  model  of  the  motor  assembly  shown  in  Figure  6 

LUMPED  NODE  MODEL  SIMULATION 

After  all  the  resistor  and  capacitor  values  of  the  lumped-node  network  model  of  the  motor  in 
Figure  6  are  known,  the  same  values  of  the  resistors  and  capacitors  can  be  used  to  simulate  the 
temperature  response  of  the  motor  parts  with  any  combination  of  heat  losses  and  boundary 
conditions.  For  this  standard  resistor-capacitor  electrical  network,  a  set  of  first-order  ordinary 
differential  equations  can  be  written  as  (Eq.  3-Eq.  13),  where  the  variable  U  represents  the  node 
temperature. 


U1-U2 

R1 


+  n^  =  /c(t) 


(4) 


U2-U3  dU 2  U1—U2 

R2  dt  R1 


-  0 


(5) 
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P4-  P8  .  P10  -  P8  PI  -  P8  .  nAU 8  ^ 

- si-+C8-«r"wr) 

(6) 

m  -b  flIQjPIO  -  P8tf9  +  ff'lQP9  +  OTfflO/lV{f> 

(7) 

PS-P8  P1S-P8  ^dUQ  „ 

- ffa — C8"ef"° 

(8) 

P13-P4  P13-P8.  _rfP18 
“RH - Rs— +a8TF",w> 

(9) 

P14-  P4  .  PIS  -  P4  .  P4-  P8  P4-  PS  ^  cfP4  „ 

Mfl  1  R: 18  1  /?8  /?4  C'f  "  ° 

(10) 

P12  -  P14  .  U7  -  P12  „„AV  12 

— rtz— +— r? — a2-ar";em 

(11) 

P9-P0  PIC  -  P9  .  _rfP9 

-RJ - RS—' tC'-3f,w 

(12) 

P6-P7  P6-P8  P9-P0  ^flPO  „ 

na  1  flop  M  *  At  ■  ° 

(13) 

P7-P12  P6-P7._dP7_ 

RJ - m-+:,-3Fm  iR7a 

(14) 

P12-P14  P14-P4  ^ 4  rfP14  ID(lA 

—m - r —*cusrm,R1* 

(15) 

This  equation  set  can  be  solved  by  standard  numerical  methods.  In  Figure  8  the  simulation 
results  of  the  lumped  node  model  are  compared  with  the  FEA  results  under  same  initial  condition 
and  boundary  conditions. 
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Figure  8.  Lumped-node  network  simulation  result 

The  results  show  excellent  accuracy  in  both  steady  state  and  transient  response  for  the  step  heat 
load  given  in  Table  3  and  Figure  4.  This  model  is  also  tested  with  pulsed  load,  where  the  thermal 
load  is  a  square  wave  form  repeated  with  time,  as  shown  in  Figure  9.  We  compared  the  resulting 
temperature  change  in  copper  winding  of  the  motor  (Figure  10). 
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Pulsed  load  definition 


Time  (s) 


0.08 

0.07 

0.06 

0.05 

0.04 

0.03 

0.02 

0.01 
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a> 


o 

Q. 


Figure  9.  Pulsed  heat  load 


Figure  10.  Winding  temperature  comparison  (FEA  and  Lumped  node) 

The  excellent  agreement  between  ANSYS  results  and  lumped-node  model  proves  that  this  model 
is  an  effective  alternative  to  computational  intensive  FEA  simulation.  The  lumped  node  model  is 
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much  faster  to  obtain  the  results  (less  than  a  minute  compared  to  20  hours  on  same  computer  for 
750  seconds  of  heat  load  profde).  It  can  also  be  easily  incorporated  into  lumped  node 
electromagnetic  motor  model  to  do  the  multi-physics  simulation. 

CONCLUSIONS 

The  lumped  node  thermal  network  is  a  feasible  approach  to  simulate  the  mission-level  EMA 
thermal  performance.  A  new  lumped  node  modeling  technique  which  is  based  on  tuning  the  R 
and  C  values  with  the  FEA  results  is  developed  in  this  work  to  overcome  the  difficulty  in 
obtaining  the  R  and  C  parameters  in  the  motor  thermal  network.  The  model  built  with  this 
method  can  be  incorporated  into  lumped  node  electromagnetic  simulation  code  to  perfonn  the 
mission-level  real-time  simulation  of  the  whole  EMA. 
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